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Exclusive decays of -B and D mesons are discussed with the focus on the theoretical 
problem how to calculate the relevant hadronic matrix elements of weak operators in 
the framework of QCD sum rules. The three lectures are devoted to 
i-O \ 1. leptonic decays and decay constants 

2. semileptonic decays and form factors 
Oh. 3. nonleptonic decays and nonfactorizable amplitudes. 

I shall introduce some of the basic concepts, describe various calculational techniques, 
and illustrate the numerical results. The latter are compared with lattice and quark 
model predictions and, where possible, with experimental data. Applications include 
the determination of V u b from B — ► nlui and B — > plvi, and an estimate of the 
phenomenological coefficient 02 for B — * J /ipK. 



1 Introduction 

Inclusive and exclusive decays of heavy flavours play a complementary role in the determina- 
tion of fundamental parameters of the electroweak standard model and in the development 
of a deeper understanding of QCD. While the theory of inclusive decays Q is well advanced, 
inclusive measurements are generally quite difficult. Conversely, exclusive decays into few- 
body final states are often easier to measure, but the theory of exclusive processes is more 
demanding and hence still underdeveloped. In view of the exciting experimental prospects 
at future bottom and charm factories, where many exclusive channels are expected to be 
measured accurately, it is pressing to make further progress on the theoretical side. 

1 See, e.g., the lectures by N.G. Uraltsev in this volume [Q. 



The difference in complexity of inclusive and exclusive observables is illuminated by 
the fact that the simple parton picture of charm and beauty decays illustrated in Tab. 1 
provides surprisingly reasonable estimates for the lifetimes and inclusive branching ratios, 
whereas no such description exists for exclusive decays. The theory of both inclusive and 
exclusive processes is based on operator product expansion (OPE) which allows to separate 
the dynamics at short and long distances. At very large scale, \i = 0(mw), charm and 
beauty decays are described by second order weak interaction involving W exchange. Since 
the momentum transfer p 2 <C one effectively has four-fermion interactions described 
by a hamiltonian of the form 

G p 

where and j' are V — A quark or lepton currents. In the nonleptonic case, evolution of 
the hamiltonian to the physical scale of the order of the heavy quark mass, /x = 0(m,Q), 
leads to considerable modifications of ( |1.1| ) due to strong interactions. These effects can be 
taken into account by OPE and renormalization group methods. The result is an effective 
hamiltonian 

Hnl = —7= Yl c i( a s> (1-2) 
V2 i 

involving a sum of local operators Oj with coefficients which can be calculated in pertur- 
bation theory as long as /i > Aq C d so that one is only dealing with interactions of quarks 
and gluons at short distances. The decay amplitudes are then given by 

A(P - X,) = %5>( a »> | Oi(//) | P) (1.3) 
V2 i 

showing the main problem, that is the calculation of the matrix elements of the weak 
operators Oi which incorporate the long-distance effects. For semileptonic decays one instead 
has to compute matrix elements of weak currents: 

A(P - W t X h ) = ^=(X h | j, | P){W l | jl | 0) . (1.4) 

This simplifies the problem. 

The usual way to proceed in inclusive decays is to assume duality of the sum over all 
possible hadronic states and the quark (and gluon) final states. The theory is further 
improved by including radiative corrections to the decay widths, quark mass effects, and 
initial bound state corrections [0. In comparison, the calculation of the hadronic matrix 
elements in ( |1.3| ) and (|1.4|) for exclusive final states Xh constitutes an ab initio nonpertur- 
bative problem. Tab. 2 illustrates typical exclusive decays considered in these lectures, and 
the corresponding matrix elements such as decay constants, form factors and nonleptonic 
two-body amplitudes. Here, the interplay between weak and strong interactions is even 
more intricate than in inclusive decays. Obviously, in order to disentangle the properties of 
weak interactions such as the charged current structure, flavour mixing and CP violation 
from exclusive processes, one must understand the impact of QCD beyond perturbation 
theory. 

Current approaches include lattice calculations, QCD sum rules, heavy quark effective 
theory (HQET), chiral perturbation theory (CHPT), and phenomenological quark models. 
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Table 1: Inclusive decays of leptons and free quarks. 
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Table 2: Exclusive decays and hadronic matrix elements. 



Decay channels 


Hadronic matrix elements 
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Each of these approaches has advantages and disadvantages. For example, quark models 
are easy to use and good for intuition. However, their relation to QCD is unclear. On 
the other hand, lattice calculations are rigorous from the point of view of QCD, but they 
suffer from lattice artifacts and uncertain extrapolations. Furthermore, effective theories are 
usually applicable only to a restricted class of problems, and sometimes require substantial 
corrections which cannot be calculated within the same framework. For example, HQET is 
very powerful in treating b — > c transitions, but a priori less suitable for b — > u transitions, 
while CHPT is designed for processes involving soft pions and kaons. 

These lectures focus on applications of QCD sum rules 0, [| |J] to exclusive B and D 
decays. Proceeding from the firm basis of QCD perturbation theory this approach sys- 
tematically incorporates nonperturbative elements of QCD. Schematically, hadronic matrix 
elements of the kind indicated in Tab. 2 are extracted from suitable correlation functions 
of quark currents at unphysical external momenta, rather than estimated directly. Thereby, 
one makes use of OPE, the analyticity principle (dispersion relations), and S-matrix uni- 
tarity. In addition, one assumes the validity of quark-hadron duality. The long-distance 
dynamics is parameterized in terms of vacuum condensates or light-cone wave functions. At 
present, these nonperturbative quantities cannot be calculated directly from QCD. There- 
fore, they are determined from experimental data. Nevertheless, because of the universal 
nature of the nonperturbative input, the sum rule approach retains its predictive power. 
Moreover, the method is rather flexible and can therefore be applied to a large variety of 
problems in hadron physics. Experience has shown that sum rules are particularly suited 
for heavy quark physics. The more recent development of sum rule techniques for exclusive 
B and D decays appears similarly promising. 

The lectures are organized as follows. The first lecture is devoted to leptonic decays and 
decay constants, the second one to semileptonic decays and form factors, and the third to 
nonleptonic decays and the problem of nonfactorizable amplitudes. I shall introduce the ba- 
sic concepts, describe various calculational techniques, and show selected numerical results, 
aiming at a compromise between the demands by theorists and the needs of experimenters. 
The results are compared with lattice and quark model predictions and, where possible, 
with experimental data. Phenomenological applications include the determination of V u t, 
from B — > irlvi and B — > plui, and an estimate of the effective coefficient a-i for B — ► J/ipK. 
In addition, I shall present a few illustrative predictions for D mesons. 

Although the charm and beauty flavours could be treated in parallel, for definiteness, I 
will usually refer to B mesons. In most cases, it is obvious how to obtain the corresponding 
results for D mesons. However, when convenient I will also use the generic notation P (V) 
for a heavy pseudoscalar (vector) meson and Q for a heavy quark. 

2 Leptonic decays and decay constants 



2 Further details can be found in a recent review by A. Khodjamirian and the present author 
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2.1 Decay width 



The leptonic B decays are induced by the weak annihilation process bu — > lv\ where I = 
e, /i, r. The relevant weak hamiltonian is given by 

H L = ^K6(^&)(T>) + h.c. (2.1) 

with T M = 7 M (1 — 75). Generically, using the decay constant fp as defined by the matrix 
element of the axial- vector current, 

(0 I qi^Q I P) = if P q» , (2.2) 
g M being the P four-momentum, one obtains the following expression for the decay width: 

T(P - in) = ^\v qQ f mpm > (1 - 4-Xfp • (2-3) 

The decay constant fp characterizes the size of the P-meson wave function at the origin, 
and therefore the annihilation probability. The proportionality to the square of the lepton 
mass mi is enforced by helicity conservation leading to a strong suppression of decays into 
light leptons. 

From the measurement of a given leptonic decay width one can directly determine the 
product of a specific CKM matrix element and decay constant. Unfortunately, because 
of multiple suppression it is difficult to measure leptonic decays. So far, only the mode 
D s — > fiis^ has been observed ||, while a bound J7J exists on the Cabibbo-suppressed decay 
D — > Jiv^. Whether or not it is possible to study B — > jiv^ and B — > fu T at future B 
factories is not completely clear ||. 

Theoretically, the decay constants fp are extremely interesting quantities. They repre- 
sent the simplest hadronic matrix elements and therefore provide crucial tests of nonper- 
turbative methods in QCD. Moreover, the decay constants are needed in sum rule analyses 
of other hadronic properties of B and D mesons, e.g., form factors. Finally, fp and fp a are 
important parameters of mixing and CP-violation in the B system. For these and other rea- 
sons which will be pointed out in the course of the lectures, accurate theoretical predictions 
on the various decay constants are very desirable. 



2.2 QCD sum rule for decay constants 

The QCD sum rule estimate of fp is based on an analysis of the two-point correlation 
function [|, g 0, [ll[ 



II(g 2 ) = i J d 4 xe iqx (0 I T{q(x)i-f 5 b(x),b(0)i~f 5 q(0)} | 0) , (2.4) 

q being the external momentum. By inserting a complete set of states with S-meson 
quantum numbers between the currents in ( |2.4j ) one obtains a formal hadronic repre- 
sentation of n(g 2 ). The term arising from the ground state S-meson is proportional to 
(0 I qi'y^b \ B){B \ bi^^q \ 0), i.e., to fp as can be seen from the relation 

m b (0 I ui<y 5 b \ B) = m 2 B f B (2.5) 



5 



which is equivalent to ( |2.2|) . The hadronic expression for ( |2.4|) can be written in the form 
of a dispersion relation: 

ntf ) = /! ^ (2.6) 



s - r 



with the spectral density 

p(s) = 5(s - ml)^& + p h (s)Q(s - s h ) . (2.7) 



m 2 b 



Obviously, the ^-function term on the r.h.s. of ( |2.7p represents the B meson, while p h (s) 
and Sq are the spectral density and threshold energy squared of the excited resonances and 
continuum states, respectively. In order to make the dispersion integral Q2.6p ultraviolet- 
finite, one actually has to subtract the first two terms of the Taylor expansion of II(g 2 ) at 
q 2 = 0. These subtraction terms are removed by Borel transformation with respect to q 2 : 



-q 2 /n=M 2 



(_ g 2)(n+l) / d 

n\ \dq 2 



B M2 U(q 2 )= lim ( — n(g 2 )^n(M 2 ). (2.8) 



With 



(^f V ' /M! ^ 



B M <-q'f = (2.10) 



and 

for k > one readily finds 

iy-y-l 4 -/*2 POO 

H(M 2 ) = H^B e -m%/M 2 + S ph( s)e s/M* ds (2 n) 

m{ hi 

We see that Borel transformation removes arbitrary polynomials in q 2 and suppresses the 
contributions from excited and continuum states exponentially relative to the ground-state 
contribution. The second point is actually the main motivation for this transformation. 

In the space-like momentum region, q 2 < 0, where one has no poles and cuts associated 
with physical states it is possible to calculate the correlation function II(g 2 ) in QCD in terms 
of quark and gluon degrees of freedom. To this end, the T-product of currents in ( [2.4|) is 
expanded in a series of regular local operators 

n( q 2 ) = Y,c d (q 2 ,p)(o \nM\o), (2.12) 

d 

d being the dimension of a given operator and fi being the renormalization scale. The 
lowest-dimensional operators (d < 6) follows: 

n d = 1, qq, G%G a ^, qa^G^q, (qT r q)(qT s q) . (2.13) 

Here, A a denotes the usual S'?7(3)-colour matrices, is the gluon field strength tensor, 
and T t stands for a given combination of Lorentz and colour matrices. While the strong 
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interaction effects at momenta larger than [i are included in the Wilson coefficients Cd{q 2 , //), 
the effects at momenta smaller than /i are absorbed into the matrix elements of the operators 
Thus, if fi ^> Aqcd the coefficients depend only on short- distance dynamics, and can 
therefore be calculated in perturbation theory, while the long-distance effects are taken into 
account by the vacuum averages (0 | fid | 0) = (fid)- These so-called condensates describe 
properties of the full QCD vacuum and are process-independent. At present, they can only 
be estimated in some crude approximations. For this reason, the vacuum condensates are 
determined empirically by fitting selected sum rules to experimental data. 

It is essential for the whole approach that at q 2 <C ml the expansion ( |2.12| ) can be cut off 
after a few terms. The reason is that the higher the dimension of fid, the more suppressed 
by inverse powers of ml — q 2 is the corresponding Wilson coefficient Cd- In leading order 
and up to d — 6, the coefficients can be calculated from the diagrams depicted in Fig. 1 as 
explained in || (see also H). In order to improve the accuracy of the OPE one can further 
include perturbative QCD corrections to the Wilson coefficients originating from hard gluon 
exchanges in the diagrams of Fig. 1. Most important are the 0(a s ) effects on the coefficient 
Co [|], ||IJ given by the two- loop diagrams of Fig. 2. 




Figure 1: Diagrams determining the Wilson coefficients in the OPE of the two-point corre- 
lation function l\2.4) : C (a), C 3 (b), C 4 (c,d), C 5 (b,e), C 6 (b,f). Solid lines denote quarks, 
dashed lines gluons, wavy lines external currents. Crosses indicate vacuum fields. 



The complete result for the Borel-transformed correlation function ( |2.4| ) obtained in JTT ] 
is given by 

s — mi)' 2 /_ 4a, „, Q .\ , s 



^,».9n 3 f 00 , (s-ml) 2 / Aa s ,. ,A 
n M 2 = — / ds K - b -L 1 + -±f(s, m 2 b ) exp 



M 2 ' 



+(- M«) + - - ^WM) 
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"i r 



b 

(a) (b) (c) 

Figure 2: Feynman diagrams of the 0(a s ) correction to the Wilson coefficient Co(q 2 ). 

l6 * a ';i-4±-^L)(«r)eM-4), (2.i4) 



where 



27M 2V AM 2 12M AJX ^' J M 2 ' 

9 „ r-Z- dt, ^ , , s 



f{s,m 2 ) = --2 f s ^ln{l-t)+\n— 2 \n 
4 jo t ml 



b s-m 2 b 

3, ml , s ml ^ s — ml ml , s 

+ -ln ^-o+ln o+— ln T ± + ~^ ln ^I ■ 2 - 15 ) 

2 s — mi s — mi s ml s — ml ml 

The 6-quark mass appearing in the perturbative contribution (first term in ( |2.14| )) is defined 
to be the pole mass, while the choice of m& in the leading-order coefficients of the contri- 
butions from the higher- dimensional operators is arbitrary. Usually, the pole mass is used 
everywhere. In order to keep the formula readable, the scale dependence of the running 
coupling and the vacuum condensates is not made explicit in the above. Furthermore, use 
has been made of the conventional parametrization for the quark-gluon condensate density, 

(qa^G^q) = m 2 (qq) , (2.16) 

and of vacuum saturation reducing four-quark condensates to squares of quark condensates 
with known coefficients [Q: 

(qT rq qT sq ) = —^-—{(TrT r )(TrT s ) - Tr(T r T s )}(qq) 2 • (2.17) 



The two representations of U(M 2 ), Q2.ll ) and (|2.14| ), yield an interesting equation which 



connects hadronic properties with QCD parameters. However, in order to determine the 
value of fs from it, one has to subtract the unknown contribution from the excited and 
continuum states to Q2.ll ). In a reasonable approximation based on quark-hadron duality, 



one may estimate the integral over the hadronic spectral density p h by the corresponding 
integral over the perturbative density. This is achieved by substituting 

p h (s)Q(s - sj}) = -lmC (s)Q(s - (2.18) 

7T 



in (|2.11|) , whereby J is treated as an effective threshold energy of the order of the mass 



of the first excited B resonance. In this approximation, continuum subtraction amounts to 
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a simple change of the upper limit of integration in ( 2.14j) from oo to Sq . The uncertainty 
from this rough procedure should not be too harmful because of the suppression of excited 
and continuum states after Borel transformation. 

Solving then the equation provided by ( [2.11| ) and ( |2.14| ) for fs, one obtains the following 
sum rule: 



p2 4 
JB m B 



3m? 



,2 c 



1 + -^-f( s i m T, 



mi 



exp 



M 2 




(2.19) 

Z,i 1V1~ \ 'ilVl' LZ.1V1 ' I I V 1V1~ I 

Since the l.h.s. of this relation is a measurable quantity, any scale dependence on the 
r.h.s. must cancel, at least in principle. In practice, this can of course only be achieved 
approximately. 



2.3 Numerical results 

The numerical input in the sum rule (|2.19 ) for fg is as follows. The pole mass of the b 
quark is taken to be 

m b = 4.7 ±0.1 GeV , (2.20) 



covering the range of estimates obtained from bottomonium sum rules while the con- 
tinuum threshold in the B channel is fixed to the range 

Sq = 35 =f 2 GeV 2 . (2.21) 

The scale fi at which a s and the condensates are to be evaluated is still ambiguous in the 
approximation considered. A reasonable choice is 

/j, 2 = 0(M 2 ) (2.22) 



characterizing the average virtuality of the quarks in the correlator (|2.4| ). Moreover, the 
products m^qq) and (^-G'J ll/ G afJ-u ) are renormalizat ion-group invariant. Using the numerical 
value of the quark condensate density at fi = 0(1 GeV) as derived from the PCAC relation 

IB 

<™ >(1 &V) = -2^k^I) * - (24 ° MCV) ' (2 ' 23) 

and the running mass rribilGeV) corresponding to the central value ( 2.20]) of the pole mass, 
one gets 

m b (qq) ~ -0.084 GeV 4 . (2.24) 

This is the estimate adopted for the present analysis together with the gluon condensate 
density determined from charmonium sum rules @: 

(y 

(—G a uu G aflu ) ~ 0.012 GeV 4 . (2.25) 

7T ^ 



The remaining parameters are assumed to be 



mg(l GeV) ~ 0.8 GeV 2 

as extracted from sum rules for light baryons fli~4 |, and 

a s (qq} 2 = 8 ■ 10~ 5 GeV 6 

as given in [flT[] . The scale dependence of m^qq) and a s (qq) 2 is negligible. The numerical 
uncertainties in the condensate estimates vary from 10 % to about 50 % or more. They 
nevertheless play only a minor role for the total theoretical uncertainty on fg. 



(2.26) 



(2.27) 



Table 3: Decay constants of B and D mesons in MeV. 



Method 


Ref. 


Sb 




Sd 


fn 3 


QCD sum rules 


this 
review 

15] a ) 


180 ± 30 
175 


210 


190 ± 20 
180 ± 10 


220 ± 10 


Lattice 


[16] 
[17] 
[18] 


175± 25 
180± 32 
1721H 


200± 25 
205± 35 
196l| 


205± 15 
221± 17 

19113 


235 ± 15 
237± 16 
2061^ 


Experiment 


1 






< 310 


241 ± 21± 30 



a) 



update of results in || [II], [L9| taking m b = 4.67 GeV, Sq = 35 GeV 2 



and m r . = 1.3 GeV, s£ = 5.5 GeV 2 . 



6 ) world average 



Next, one has to determine the range of values of the Borel parameter M 2 for which the 
sum rule ( |2.19|) can be trusted. On the one hand, M 2 has to be small enough such that 
the contributions from excited and continuum states are exponentially damped and the 
approximation Q2.18 ) is acceptable. On the other hand, the scale M 2 must be large enough 
such that higher-dimensional operators are suppressed and the OPE converges sufficiently 
fast. Furthermore, the sum rule should be stable under variations of M 2 in the allowed 
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window. Whether or not these requirements can be met has to be investigated for each sum 
rule separately. In the case at hand, one finds that at 



3 GeV 2 < M < 5.5 GeV 2 (2.28) 

the necessary conditions are fulfilled: the quark-gluon condensate contributes less than 15 
%, and the four-quark operators a negligible amount, while the excited and continuum 
states contribute less than 30 %. Moreover, changes insignificantly when M 2 varies in 
the range (ggp . 

The predictions on fs and fs a derived from the sum rule ( |2.19| ) with the input specified 
above are summarized in Tab. 3 together with the corresponding results for D mesons 
obtained with m c = 1.3 ± 0.1 GeV, s» = 6 T 1 GeV 2 , and 1 GeV 2 < M 2 < 2 GeV 2 . 
The theoretical uncertainties are dominated by the uncertainties in the quark masses and 
continuum thresholds. For instance, the ranges of values of fs and fo expressed by the 
errors in the first row of Tab. 3 have been determined by varying uiq and Sq simultaneously 
such that the stability of the sum rule is preserved. The correlation is indicated by the 



reversed ± sign in ( 2.21 ). It should also be emphasized that the 0(a s ) corrections have a 



sizeable effect. Without these corrections the decay constants would be considerably smaller: 

J B = f B (a s = 0) = 140 ± 30 MeV , (2.29) 

J D ee f D {a s = 0) = 170 ± 20 MeV . (2.30) 

For comparison, Tab. 3 also shows recent lattice results and the available experimental 
data. The mutual agreement within the uncertainties of the two theoretical approaches is 
satisfactory. On the other hand, the data are just beginning to challenge theory. 



3 Semileptonic decays and form factors 
3.1 Distribution in momentum transfer 

The weak hamiltonian (|2.1|) also applies to the semileptonic decay B — > ixlvi. The appro- 
priate hadronic matrix element is parametrized by two independent form factors: 



(TT(q)\u^b\B(p + q)) = 2f + {p 2 )q, + (f+(p 2 ) + r{p 2 )) Pfl 



(3.1) 



with p + q, q and p being the B and ir four-momenta, and the momentum transfer, respec- 
tively. Frequently, in particular in nonleptonic decays, one also uses the scalar form factor 
f° defined by 



/V) = f + (p 2 ) + 



p 



ni 



mi 



-r( P 2 



(3.2) 



The distribution of the momentum transfer squared in B — > irlvi is given by 



dT_ = G 2 \V ub \ 2 (p 2 -mjf^JEl-m 2 
dp 2 ~ 24tt 3 



m 



1 + ^)^(^-0 f + (p* 



+ 



p 4 m 2 B 
3m 2 , 



B 



,2\2 



(3.3) 
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where E n = {m? B + m\ — p 2 )/2m_g is the pion energy in the B rest frame. 

For I = e, /i, the form factor f° plays a negligible role because of the smallness of 
the lepton masses. These channels can thus be used to directly determine the product 
\V u bf + (p 2 )\, while the heavy mode B — > 7ifu T offers a possibility to test the scalar form factor 
f°{p 2 )- Recently, the decays B — > 7r/z/; and -B — > pZz/j have been observed by CLEO |20 



If the relevant form factors can be calculated with sufficient accuracy these measurements 
eventually provide interesting alternatives to the determination of V u b from inclusive b — > u 
transitions. In these lectures, the sum rule method is advocated. It is one of the virtues of 
the sum rule calculations, in comparison to other approaches, that they can be tested and 
consistently improved by investigating the analogous D-meson form factors. The already 



existing experimental data on D — > hlui, h = K, K*,ir, p,rj and I = e,/i |21| have not yet 
been fully exploited in this respect. 

3.2 Light-cone vs. short-distance expansion 

For processes involving a light meson, e.g., a n, K, or p there is an interesting alternative 
to the short-distance OPE of vacuum- vacuum correlation functions in terms of condensates, 
namely the expansion of vacuum-meson correlators near the light-cone in terms of meson 
wave functions p2[ ^3], [24) . The latter are defined by matrix elements of composite operators 
and classified by twist, similarly as the deep- inelastic structure functions. The light-cone 



variant of QCD sum rules has been suggested in |25|, ^6], pvjj . In this approach, the light 
mesons are described by means of light-cone wave functions, whereas the heavy meson 
channels are treated in the usual way: choice of a generating current, dispersion relation, 
Borel transformation and continuum subtraction. 

The general idea is explained in the following using the B — > 7r matrix element ( |3.1| ) as 
an example. The detailed derivation of the light-cone sum rules for / + and /~ is reported in 
and reviewed in ||. The central object is the vacuum-pion correlation function 

F,(p,q) = % [ d 4 xe tpx (7r(q) | T{u(x)^b(x), b{0)i l5 d(0)} | 0) (3.4) 



= F(p 2 , {p + qf)q^ + F(p 2 , (p + qf)p^ . 

From the comparison with (|3.1| ) it is clear that the two invariant functions F and F have 
some kind of relation with the form factors / + and f + + f~ , respectively. Since the pion is 
on-shell, q 2 = m 2 vanishes in the chiral limit adopted throughout this discussion. For the 
momenta in the bd and ub channels one requires (p + q) 2 ml and p 2 < m 2 — 2mbXi X being 
a mfe-independent scale of order Aqcd- This assures that the b quark is far off mass-shell. 
Contracting the 6-quark fields in (^.4[) and inserting the free 6-quark propagator, one gets 



+ k»(7r(q)\u(x) ltilul5 d(0)\0)) . (3.5) 

This contribution is diagrammaticly depicted in Fig. 3a. It is important to note that the 
path-ordered gluon operator 

Pexp ^ig s J da x fl A tJ '(ax)'^ (3.6) 
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necessary for gauge invariance of the matrix elements in ( |3.5|) , is unity in the light-cone 
gauge, x^A^ = 0, assumed here. 




Figure 3: Diagrammatic representation of the correlation function (oJ^j: terms involving 
two-particle (a) and three-particle (b) wave functions. Solid lines represent quarks, dashed 
lines gluons, wavy lines external currents, and ovals light-cone wave functions of the pion. 

Let us first consider the matrix element 

(7r(g)|iZ(ar)7 M 7 5 d(0)|0), (3.7) 

and expand the bilocal quark-antiquark operator around x = 0: 

u(x) W d(0) = £ - } u(0)(D -x) r 1,1^(0) . (3.8) 

r T '• 

The matrix elements of the local operators can be written in the form 

(n(q)\u D ai D a2 ... D ar 7^7s^|0) = (i) r ' q^q ai q a2 ---q ar M r + ... , (3.9) 

where the ellipses stand for additional terms containing the metric tensor g ai a k in all possible 
combinations. Note that the lowest twist f\ in (|3.9| ) is equal to 2. Substituting (|3.8| ) and 
( |3.9| ) in the first term of ( |3.5| ), integrating over x and k, and comparing the result with fl3.4|) , 
one obtains 

F{p\ (p + qf) = t^- 2 £ CM r (3.10) 

with 

{= jfe^) = (? + «)' -/ . (3.11) 
ml — p z ml — p z 



3 Twist is defined as the difference between the canonicai dimension and the spin of a traceless and 
totally symmetric local operator. 
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Obviously, if the ratio £ is finite one must keep an infinite series of matrix elements in ( |3.10 ). 
All of them give contributions of the same order in the heavy quark propagator I /(ml — p 2 ), 
differing only by powers of the dimensionless parameter £. In other words, the expansion 
( |3.8| ) is useful only in the case £ — > 0, i.e., p 2 — > (p + q) 2 or equivalently, when the pion 
momentum vanishes. In this case, the series in ( |3.10|) can be truncated after a few terms 
involving only a manageable number of unknown matrix elements M r . However, generally, 
when p 2 ^ (p + q) 2 one has to sum up the infinite series of matrix elements of local operators 
in some way. 

One possible solution of this problem is provided by the light-cone expansion of the 
bilocal operator (|3.8|) . For the matrix element ( |3.7|) the first term of this expansion is given 
by 

(Tr(q)\u(x) llll5 d(0) \0) x i =0 = -iq^U t due^^u) . (3.12) 

Jo 

The function ip n (u) is known as the twist 2 light-cone wave function of the pion. It represents 
the distribution in the fraction of the light-cone momentum q + q 3 of the pion carried by 
a constituent quark. As can be seen by putting x = in the above, <p n (u) is normalized to 
unity. Substitution of Q3.12j ) in (|3.5|) and integration over x and k yields 

F(p\ (P + qf) = m b U f 1 . (3.13) 

Jo mf — (p + uq) z 

We see that the infinite series of matrix elements of local operators encountered before in 
( |3.10| ) is effectively replaced by a wave function. Formal expansion of ( |3.13| ) in q, 

00 (2r> ■ nY r 1 
F(p\ (P + qf) = mU E A P / du vT^{u) , (3.14) 

r=0 ' h " ' 



and comparison with (|3.10|) yields the relation between the matrix elements M r defined in 
(|3.9|) and the moments of (p n (u): 

M r = —if-K / du litpJu). (3.15) 
Jo 

Including the next-to- leading terms in x 2 , the light-cone expansion of the matrix element 
is given by 



(7r(g)|n(a;)7 M 75d(0)|0> = -iq^f n / due*"** (cp w (u) + x 2 9l (u) 

+U U» - — J C due iu ^g 2 (u), (3.16) 
\ qx ) Jo 

where g± and g 2 are twist 4 wave functions. Proceeding to the second term in (|3.5|) and 
using the relation 

7/^75 = 9^,15 ~ i<VTk , (3-17) 
one encounters two further matrix elements: 

(ir{q) | u(x)i l5 d(0) | 0) = f nf ^ [ du e iuqx (p p (u) (3.18) 
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and 



(ir(q) | u(x)a lxu j 5 d(0) | 0) = i{q^x u - q u x^ 



du e iuqx Va (u) 



(3.19) 



6 Jo 

with p n = m 2 /(m u + m^). Only the leading terms of the expansion are considered here. 
They have twist 3 and are parameterized by the wave functions (p p and (p a . 

Beyond twist 2, one should also take into account the higher-order term illustrated 
graphically in Fig. 3b, where a gluon is emitted from the heavy quark line. This correction 
is determined by quark-antiquark-gluon wave functions |31|. Explicit expressions for this 
contribution can be found in (gBJ, 29| and also in ||. Below, this term is denoted by Fq- 
Gluon radiation from the light quark lines in Fig. 3a is effectively included in the twist 3 
and 4 wave functions as was shown in [^, |32| \ 



. Components of the pion wave function 
with two extra gluons, or with an additional qq pair are neglected. Including all operators 



up to twist 4 one obtains the following results for the two invariant functions defined in 



± 

F(p 2 ,(p + q) 2 ) = Uj 



F(p 2 ,(p + qf 



du 



mi 



(p + uq) 2 



m b ip w (u) 



uip p (u) + - \ 2 + 

2ug 2 {u) 

m l ~ (P + U( l) 2 
+ F G (p 2 ,(p + q) 2 ) 
1 

du 



+ Hit 



+ m b 



p 2 + m 2 



mi 



(p + uq) 2 



ip a {U) 



imi 



m\ — (p + uq) 2 ) 2 



9iW 



dvg 2 (v) 



(3.20) 



mi 



(p + uq) 2 



mi 



mi 



(p + uq) 2 



+ 



6u 
2m b g 2 {u) 



ml — (p + uqy 



(3.21) 



We see that every power of x 2 in the light-cone expansion of the integrand in ( |3.4| ) leads to 
an additional power of ml — (p + uq) 2 in the denominators of the coefficients. This justifies 
the neglect of higher-twist operators, provided (p + q) 2 and p 2 are sufficiently smaller than 
m\. It is also interesting to note that there are no contributions from twist 2 and from 
three-particle wave functions to F. 



3.3 Light-cone wave functions 

Before proceeding with the derivation of the light-cone sum rules for the form factors / ± , it 
may be useful to make a few more remarks on the light-cone wave functions. Firstly, these 
functions are universal, a property which is essential for the whole approach, similarly as the 
universality of the vacuum condensates in the case of short-distance sum rules. Secondly, 
the asymptotic form of the wave functions is given by conformal symmetry and asymptotic 
freedom of QCD [^4], |33|], while the nonasymptotic features incorporate long-distance quark- 
gluon interactions. Thirdly, the factorization of collinear logarithms generated by gluon 
radiation and the renormalization procedure induce scale dependences. Usually, a common 
scale fj, is chosen for simplicity. 
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In leading-order approximation (LO), the twist 2 wave function ip n (u, fi) obeys the 
Brodsky-Lepage evolution equation [f22|j: 



d In fi 



dwV(u, w)ip n (w, /i) 



where 



V(u, w) 



a s (n)C f 



TV 



1 - U 

1 — w 



1 + 



1 



u — w J 



u-w) + — [l + 
w 



w — u) 



) B(w-u) 



and the + regularization is defined by 



(3.22) 



(3.23) 



R(u, w) 



R(u, w) — 5(u — w) R(v, w)dv 



(3.24) 



The evolution equation effectively sums up the leading logarithms to all orders. It can be 
solved by expanding <p n (u, [/,) in terms of Gegenbauer polynomials 



(f n (u, fi) = 6u(l -u) 1 + a 2i(/i)C , 2/ 2 (2w - 1 



(3.25) 



in which case the coefficients ci2i(/f) are multiplicatively renormalizable [22]. The scale 
dependence is given by 



a n (jd) = a n (n ) 



' OLs(fJ>) 



In /b 



(3.26) 



where b = 11 — 2nj/3 is the 1-loop coefficient of the QCD beta function, rif being the 
number of active flavours, and 

7n = C F 



(n + l)(n + 2) 



+ 4 E 



\k=l 



k 



(3.27) 



are the anomalous dimensions |p4] . We see that a n (/i), n > 2 vanishes for /i — > oo. There- 
fore, these terms describe the nonasymptotic effects in the wave function ( |3.25| ). The /i- 
independent term is the asymptotic wave function, ip n (u, /j, = oo) = 6u(l — u). As already 
pointed out, the normalization is such that 



duipniu, ji) 



(3.28) 



The initial values of the nonasymptotic coefficients can be estimated from two-point 
sum rules [^4| for the moments J u n ip 7T (u, fi)du at low n. The nonperturbative information 
encoded in the quark and gluon condensates is thereby transmuted into the long-distance 
properties of the wave function. Alternatively, one can determine the coefficients directly 
from light-cone sum rules for known hadronic quantities such as the irNN and uipir couplings. 
Here, the estimates at fi = 0.5 GeV given in |26| are adopted: 



c^Guo) = g, «4(a*o) = 0.43 . 



(3.29) 
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The corresponding Gegenbauer polynomials are given by 



cT(2u-l) 



Q 3/2 (2m-1) = — [21(2w-l) 4 -14(2«-l) 2 + l] . 



2 
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[5(2u - l) 2 - 1] 



(3.30) 



On the basis of the approximate conformal symmetry of QCD it has been shown |33| that 
the expansion ( |3.25| ) converges sufficiently fast so that the terms with n > 4 are negligible. 

For completeness, the twist 3 and 4 two-particle light-cone wave functions of the pion 
are listed below B6jl: 



1 + B 2 (ji) l(3(u - uf - 1) + B 4 (^(35(u - uf 

Z o 



30(w-w) 2 + 3) 



6mm 



1 + C 2 (w)-(5(m -u 
15 



r.\2 



+ c 4 (fi)—(2i( u - uf - 14(« - uf + r 



5 1 6 

gx{u,ii) = -5 2 (/j,)u 2 u 2 + -e(fx)8 2 (fx)[uu(2 + 13uu) + 10u 3 \nu(2-3u + -u 

6 



+ 10w 3 lnw(2-3u + -w 2 )] 

5 



10 



5 2 (fi)uu(u — u) 



2„.2 



-5 (n)u u 



(3.31) 

(3.32) 
2 ) 

(3.33) 
(3.34) 

(3.35) 



with u = 1 — u and 



#2 = 30^-, B± = ^^^-(4a; 2i o - - 2u;i )0 ) 



fJ'irfir 



2 ^Trfn 



(■■ 



-^(5-^1,0), C 4 = ^-— ^-(4^2,o -wi,i) . 

«7T JtT 2 10 //„./„. 



The additional parameters appearing in the above are numerically given by 

iW(lGeV) = 0.0035 GeV 2 , 

w 1)0 (lGeV) = -2.88 , cj 2i0 (lGeV) = 10.5 , ^(lGeV) = , 
5 2 (lGeV) = 0.2 GeV 2 , e(lGeV) = 0.5 . 
For the quark- ant iquark-gluon wave functions the reader is referred to and 



(3.36) 



(3.37) 



3.4 Sum rules for heavy-to-light form factors 

Now we are ready to proceed with the discussion of the B — > 7r transition form factors. In 
order to determine f ± {p 2 ) from the correlation function F^(p,q) derived in (|3.4j) , we apply 
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QCD sum rule methods to the 5-meson channel following essentially the same steps as in 
the derivation of the sum rule for f B in sect. 2.2. The hadronic representation of fl3.4| ) is 
obtained by inserting a complete set of intermediate states with I?-meson quantum numbers 
between the currents on the r.h.s. of ( |3.4| ). Using the matrix elements ( |2.5| ) and ( |3.1| ), and 
representing the sum over excited and continuum states by a dispersion integral with the 
spectral density p h = p h q fl + p h Pfj,, one obtains the following relations for the invariant 
amplitudes F and F: 



2m 2 B f B f+(p 2 ) 
m b (m B — (p + qf 



+ 



p h (p 2 , s)ds 
s — (p + q) 2 



(3.38) 



jV, (p + ,).) = TiMm±rm + J ds . w. •) 



m b (m 2 B — (p+ qf 



s - (p + q) 2 



(3.39) 



Similarly as in ( 2.11 ), the integrals over p h and p h are approximated by corresponding 
integrals over the spectral densities calculated from the light-cone expansion of ( |3.4j ) using 
the substitution 

1 



—lmF(p , s)0(s — s 



71 



) ' 



(3.40) 



and the analogous relation for ph and ImF. The calculation of the imaginary parts of F 
and F from the functions given in ( |3.20|) and ( |3.21|) , respectively, is outlined, e.g., in 0. 
With the above duality approximation it is again rather straightforward to subtract the 
continuum contributions from the l.h.s. of ( |3.38| ) and ( |3.39| ). After performing the Borel 
transformation in (p + q) 2 , one arrives at the sum rules [28, |2^ 



fsf + {p 2 



fnmi I III 
7T exp 



2m 



B 



X ip n (u) + 




Am 2 gi {u) 2 



uM 2 



g 2 {v)dv 1 + 



m 2 + p 2 



u 



M 2 



+ t + {s»,p 2 ,M 



mp 2 ,M 2 ) + ^S+(p 2 ,M 2 ) 



(3.41) 



and 



fs{f + {p 2 ) + /" (P 2 )) - f ^ mb — (mB 



ml GXP \M 2 





m 2 — p 2 (l 



uM 2 



+ t ± (s*,p 2 ,M 2 



p n uM 2 



u 



(3.42) 
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The sum rule for the scalar form factor f° follows from the relation ( |3.2|) . In the above, t + 
and are numerically small corrections left from continuum subtraction, and f@ denotes the 
contribution from the three-particle wave functions. The explicit expressions of these terms 
are given in [^, |30|. Also indicated is the perturbative QCD correction 5 + to the leading 
twist 2 term in / + which has been calculated only recently [|34[ 35]. Other perturbative 
corrections are still unknown. 

Concerning the 0(a s ) correction to the leading twist 2 result, a few comments are in 
order. Hard gluon exchanges give rise to radiative effects in the correlation function fl3.4|) 
which can be calculated perturbatively. The first-order diagrams are shown in Fig. 4. The 
expression ( |3.13|) of the invariant function F can be more generally written as a convolution 



of a hard scattering amplitude T with the twist 2 wave function: 

F(p 2 , (p + q) 2 ) = U j dupn(u, /i)T(p 2 , (p + q) 2 , u, fi) , (3.43) 
J o 

where T is determined by a perturbative series 

T(p 2 , (p + q) 2 , u, /,) = T (p 2 , (p + q) 2 , u) + ^Up 2 , (p + q) 2 , u, /*) + 0(a 2 s ) . (3.44) 

The LO approximation is given by the zero-order amplitude 

Up 2 , (p + qf, u) = — m \ (3.45) 

mi — p z (l — u) — (p + gj^-u 

convoluted with the wave function obeying the LO Brodsky-Lepage evolution equation 
( 3.22|) . In next-to- leading order (NLO), the hard amplitude Q3-44|) includes the first-order 



term T 1; while the wave function (f^u, n) solves the NLO evolution equation [Kffl. For cal- 



culational details and analytic results one may consult p4| , |35| . To date, the NLO program 
is completed only for the twist 2 contribution to the invariant amplitude F. 

It is important to note that the correlation function (|3.4j) actually leads to sum rules for 



products of decay constants and form factors as made explicit in (|3.41| ) and ( 3.42|) . When 



dividing out the decay constants, consistency requires to take both the sum rule and fp in 
LO or in NLO. This appears obvious, but has not always been respected in the literature. 
Some of the discrepancies in the calculations originate just from mixing orders. In fact, as 
will be quantified later, the remaining radiative corrections to f + {p 2 ) turn out to be rather 
small because of a cancellation of the corrections to the sum rules ( |3.41| ) for /s/ + (p 2 ) and 



(TO for f B . 

As a final remark, the light-cone sum rules provide a unique possibility to investigate 
the heavy-mass dependence of the heavy-to-light form factors. To this end, one employs the 
following scaling relations for mass parameters and decay constants: 

m P — m Q + A , Sq =m 2 Q + 2m Q uj , M 2 = 2m Q r , (3.46) 

/„ = fp/Vrnj, (3-47) 

where in the heavy quark limit A, ljq, t, fp are mQ-independent quantities. With these 
substitutions it is straightforward to expand the sum rules (|3.41|) and ([3.42 ) in vtlq. In both 



cases, the light-cone expansion in terms of wave functions with increasing twist is consistent 
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Figure 4: Diagrams of the 0(a s ) corrections to the correlation function ($.4 )- 

with the heavy mass expansion, that is the higher-twist contributions either scale with the 
same power of mq as the leading- twist term, or they are suppressed by extra powers of tuq. 

The asymptotic scaling behaviour of the form factors turns out to differ sharply at small 
and large momentum transfer |3(| . At p 2 = 

/ + (0) = /°(0)~m Q 3/2 , (3.48) 

/ + (0) + r(0)~m Q 3/2 , (3.49) 
whereas at p 2 = rriq — 2m,QX, X being independent of itiq, 

f + (p 2 )~m¥ , (3-50) 
f + (p 2 ) + f-(p 2 )^f°(p 2 )~m Q 1/2 . (3.51) 

The sum rules thus nicely reproduce the asymptotic dependence of the form factors on the 
heavy quark mass derived in |S7| , |3S|1 for small pion momentum in the rest frame of the P 
meson. In addition, the sum rules allow to investigate the opposite region of large pion 

4 see also ref. 
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momentum where neither HQET nor the single-pole model can be trusted. Clearly, from 
the sum rule point of view it is expected that the excited and continuum states become 
more and more important as p 2 — > 0. This is reflected in the change of the asymptotic mass 
dependence. Claims in the literature which differ from fl3.48| ) to ( [3.51]) are often based on 



the pole model and therefore not reliable. A similar analysis has been carried out in |39| for 
B — > p form factors with essentially the same conclusions. 



3.5 Numerical analysis 

For most of the following illustrations, it suffices to work in LO. Strictly speaking, since 
the NLO order corrections to + /~) are unknown, it is also more consistent to work 

in LO when combining or comparing / + and f + + f~ ■ Therefore, if not stated otherwise, 



the values of m&, Sq, and /b(q; s = 0) specified in fl2.20| ), fl2.21Q and ( |2.29| ), respectively, are 



used. The pion wave functions are as given in sect. 3.3 with the coefficients scaled to 



fi b = ym^-mg ~ 2.4 GeV . (3.52) 

With this input we have checked that for M 2 = 10 ±2 GeV 2 the twist 4 corrections are very 
small, and that the continuum contributions do not exceed 30%. Up to p 2 ~ 18 GeV 2 , the 
sum rules are also quite stable under variation of M 2 . At larger momentum transfer the sum 
rules become unstable, and the twist 4 contribution grows rapidly. This is not surprising 
because the light-cone expansion and the sum rule method are expected to break down as 
p 2 approaches m 2 . 




p 2 [GeV 2 ] 



Figure 5: B — > 7r form factors obtained from light-cone sum rules in LO. 

The momentum dependence of the form factors / + and / + + /~ can be seen in Fig. 5. 
Specifically, at zero momentum transfer we predict 

/+(0) = 0.30, 
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/ + (0) + r(0) = 0.06. (3.53) 

If the known 0(a s ) corrections are included in the sum rules for /s/ + (p 2 ) and one gets 
instead 

/ + (0) = 0.27. (3.54) 



f + (p 2 ) 






M 2 [GeV 



12 4.6 4.7 4.8 33 35 37 



m b [GeV] 



Figure 6: Variation of the sum rule prediction on f + {p 2 ) with the Borel parameter (a), the 
b-quark mass (b), and the continuum threshold Sq (c). Considered are three typical values 
of momentum transfer: p 2 = (solid), p 2 = 8 GeV 2 (long -dashed), and p 2 = 16 GeV 2 
(short-dashed). 

Allowing the parameters to vary within the ranges specified one can estimate the uncer- 
tainties in the sum rule predictions. The main sources of uncertainty are discussed below: 

(a) Borel mass parameter 

The dependence of / + on M 2 is illustrated in Fig. 6a. In the allowed interval of M 2 , f + 
deviates by only ±(3 to 5) % from the nominal value, depending slightly on p 2 . 

(b) 6-quark mass and continuum threshold 
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Fig. 6b shows the variation of / + with mj keeping all other parameters except fixed. 
The value of is taken from the sum rule ( |2.19| ) after dropping the 0(a s ) corrections. The 
analogous test for s$ is performed in Fig. 6c. If Sq and m& are varied simultaneously such 
that one achieves maximum stability of the sum rule for the change in / + is negligible 
at small p 2 and only about ± 3 % at large p 2 . 

(c) higher-twist contributions 

No reliable estimates exist for wave functions beyond twist 4. Therefore, we use the magni- 
tude of the twist 4 contribution to / + as an indicator for the uncertainty due to the neglect 
of higher-twist terms. From Fig. 7a we see that the impact of the twist 4 components is 
comfortably small, less than 2 % at low p 2 and about 5 % at large p 2 . This suggests a 
conservative estimate of ± 5 % uncertainty due to unknown higher-twist. 

(d) light-cone wave functions 

In order to clarify the sensitivity of / + to the nonasymptotic effects in the wave functions, 
the corresponding coefficients are put to zero and the result is compared, in Fig. 7b, to the 
nominal prediction. The change amounts to about -10 % at small p 2 and +10 % at large p 2 , 
while the intermediate region around p 2 ~ 10 GeV 2 is very little affected. Since this exercise 
is rather extreme, the real uncertainty from the errors in the nonasymptotic coefficients is 
considerably smaller. 

(e) perturbative corrections 

The 0(a s ) corrections to the leading twist 2 term in the sum rule Q3.41 ) for fsf + and to 



the sum rule Q2.19] ) for both being about 30 %, cancel in the ratio. The net effect is the 



surprisingly small correction to / + shown in Fig. 8. This result [3^, |35| eliminates one main 
uncertainty. Unfortunately, the perturbative corrections to the higher-twist terms are not 
known. Considering that the twist 3 terms contribute about 50 % to the sum rule for 
and assuming again radiative corrections of about 30 % but no cancellation, the remaining 
uncertainty is 15 %. 

In summary, the present total uncertainty in the normalization of f + {p 2 ) is estimated 
to be about 25 %, if the individual contributions (a) to (e) with the exception of (d) are 
added up linearly, and about 17 % if they are added in quadrature as is often done in the 
literature. Ultimately, this may be reduced to about 15 to 10 %. In addition, there is a 
shape-dependent uncertainty from (d) of less than -10 % at low and +10 % at high p 2 . 
Moreover, in the integrated width the latter uncertainty averages out almost completely. 
The uncertainties from (a) to (d) on / + + /~ are of comparable size, while the effect of 
radiative corrections is still unknown. 

Finally, it is interesting to compare the predictions on the B — > tt form factors from the 
light-cone sum rules with the results of other approaches. Such a comparison is presented in 
Fig. 9 and 10. Although there is rough agreement on the normalization at low momentum 
transfer, the light-cone sum rules predict the fastest increase of / + with p 2 . 

The sum rules (|3.41|) and ( p.42|) for the B — > tt form factors are formally converted into 
sum rules for the D — > tt form factors by replacing b with c and B with D. The input 
parameters are taken over from the calculation of fu in sect. 2.2. In addition, one has to 
rescale the wave functions from //& ~ 2.4 GeV to /i c ~ 1.3 GeV. The allowed Borel mass 
window is 3 GeV 2 < M 2 < 5 GeV 2 . With this choice, the D — > tt form factor at p 2 = is 
predicted to be |29|, |30 | 



/+(0) = 0.68. (3.55) 
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Figure 7: B — > 7r /orm factor f + : (a) individual contributions of twist 2 (dashed), 3 (dash- 
dotted), and 4 (dotted), and the total sum (solid); (b) wave functions with (solid) and without 
(dashed) nonasymptotic corrections. 
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4 8 12 16 



p 2 [GeV 2 } 

Figure 8: B — ► 7r /orm factor f + (p 2 ) in leading twist 2 approximation: LO (dashed) in 
comparison to NLO (solid). 



0(a s ) corrections are not included here. The momentum dependence of / + in the range 
0(1 GeV 2 ) is shown in the next section, together with an extrapolation to higher 



p 2 < m 2 c 
p 2 . 



Other important applications of light-cone sum rules include the estimate of the B — ► K 

J/ipK amplitude, the 
K*~f g3 



form factor 28 



which determines the factorizable part of the B 
prediction of the matrix element of the electromagnetic penguin operator for B 
and the calculation of the B —>■ p form factors f3*5fl . In the latter work, essentially the same 
procedure is applied as outlined above. However, the relevant weak currents and light-cone 
wave functions are very different. 



3.6 Single pole approximation and the B*Btt and D*Dtt couplings 

As repeatedly stressed, the light-cone sum rules derived in sect. 3.2 and 3.4 are only valid 
in the momentum range p 2 < rriq — 2itlqX- In order to calculate the integrated semileptonic 
width T(P — > hlv\) one must therefore find another way to estimate the form factors in the 
high momentum range m 2 Q — Ivtlqx < p 2 < (mp — m^) 2 . Near the pole p 2 = m 2 / of the 
vector ground state V it is justifiable to use the single pole approximation 

Here, fy is the decay constant of the vector state defined by the matrix element 

(0 | qi^Q | V) = myfve, , (3.57) 
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and gvp-K is the strong coupling between the heavy vector and pseudoscalar mesons and the 
pion defined by 

(Vn | P) = -g VP M ■ e) , (3.58) 

e M being the V polarization vector. 

The B*Bn and D*Dtt couplings have been studied with different variants of sum rules 



and in a variety of quark models. In p9| , a light-cone sum rule for gB*B-K has been suggested 
which is derived from the same correlation function ( ft.4|) as the B — > tt form factors. 
Hence, the nonperturbative input is exactly the same in both calculations . The key idea is 
to write a double dispersion integral for the invariant function F(p 2 , (p + q) 2 ). Inserting in 
complete sets of intermediate hadronic states carrying B and B* quantum numbers, 



respectively, and using the matrix elements ( |2.5|) , ( |3.57|) and (|3.58| ) one obtains 



*V. (? + «)') = . f^/Tv' rr + / 7 i' h %f As ; Ag ; (3.59) 

m b {p z - m B ,)((p + qy - mf B ) Jn \S\ - p z ){s 2 — (p+ q) ) 
+ {subtractions) . 

The first term is obviously the ground-state contribution containing the B*Bir coupling, 
while the spectral function p h (si, S2) represents higher resonances and continuum states in 
the B* and B channels. The integration region in the (si, S2) - plane is denoted by S. The 
subtraction terms are polynomials in p 2 and/or (p+q) 2 which vanish by Borel transformation 
of ( |3.59|) with respect to both variables p 2 and (p+q) 2 - The resulting hadronic representation 



of F is thus given by 

F(M 1 ,M 2 ) = B M 2B M 2F(p ,{p + q) ) = exp 

1 2 m b 



Ml M'. t 
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Figure 10: The B — ► 7r /orm factor f°: light-cone sum rule prediction (solid) in comparison 
to the quark model result from $J\1 (dashed) and the sum rule estimate from (dash- 
dotted). 



+ 



exp 



si 



Mi 



^2 

Ml 



s 2 )dsids 2 



(3.60) 



where M 2 and M| are the Borel parameters associated with p 2 and (p + g) 2 , respectively. 
The same transformation is also applied to F(p 2 , (p + g) 2 ) given in fl3.20|) : 



F(M 2 ,M 2 ) = m b f n tp 7r (u )M 2 e~^ + (/^her twist) 



(3.61) 



with u = Ml/ '(Ml + Mf) and M 2 = Af?Mf/(Af? + Mf). The derivation of the complete 
result and, in particular, the subtraction of the continuum is quite tricky. Further explana- 
tions can be found in |J. It turns out that the subtraction greatly simplifies if M 2 = Mf , 
i.e., Mo — 1/2- With this choice, the following sum rule is obtained from ( |3.60|) and (|3.61|) 
after continuum subtraction: 



f B f B* 9 B* B 7T 



m B rriB* 







|m 2 


e "m 1 — e m 1 



</v(wo) + — (uoip p (u ) + ~<At(wo) + \ u v-P-{ u q) 
nit, V 3 6 du 



2/3tt r G/ \ 



+ e 



^7716 



y? CT (w ) + 2u g 2 (u ) 



Ami 
If 2 



(3i(«o) 



G 2 (u )) + I?\ 



u , 



(3.62) 



Here. 



uo=l/2 

and ifiua) are contributions from the three-particle amplitude shown in Fig. 



3b and given formally in |29|, Since G-parity implies g 2 (uo) 



at m = 1/2, 



these terms can be dropped in the above sum rule due to the choice M{ = M£ 
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One of the numerically crucial parameters is the leading twist 2 wave function ip n at the 
symmetry point uq = 1/2. The same parameter also enters the sum rules for other hadronic 
couplings involving the pion. For the following predictions the value 



<p v (u = -) = 1.2 ± 0.2 (3.63) 



obtained from the light-cone sum rule for the pion-nucleon coupling |26] has been used. This 
value is consistent with the coefficients 612,4 given in ( |3.29| ). The values of the remaining 
parameters are as in the calculation of the form factor / + in sect. 3.5. With this input the 
sum rule ( |3.62| ) yields 

f B fB*g B *Bn = 0.64 ± 0.06 GeV 2 . (3.64) 
Furthermore, using the sum rule estimate (p. 29 ) for f B , and 

f B * = f B * (a s = 0) = 160 ± 30 MeV (3.65) 

derived from a similar two-point sum rule for f B * , one gets 

g B . Bw = 29 ± 3 . (3.66) 

Note that the couplings of the different charge states are related by isospin symmetry: 

9b*b-k = 9b*o b -tt+ = -V^b-obOttO = V2g B »- B - 7T om = —g B *-s°ir- ■ (3.67) 

The uncertainties quoted above reflect the variation of the results with the Borel mass in 
the allowed window 6 < M 2 < 12 GeV 2 . In this window the excited and continuum states 
contribute less than 30 % and the twist 4 corrections do not exceed 10 %. There are other 
sources of uncertainties. In particular, the above predictions can and should be improved 
by calculating radiative gluon corrections. 

The sum rule ( |3.62| ) for g B * b-k is translated into a sum rule for go'Dir by formally changing 
b to c, B to D, and B* to D* . With (|3.63|) and the input parameters from the calculation 
of the D — > 7r form factor, one finds 

f D fD*9D*D-x = 0.51 ± 0.05 GeV 2 . (3.68) 

Use of ( |2.30|) for f D and of the analogous sum rule estimate 

fjj. = fjj. (a s = 0) = 240 ± 20 MeV (3.69) 

yields 

9d*Dtt = 12.5 i 1.0 . (3.70) 

The variation with M 2 in the allowed interval, here, 2 < M 2 < 4 GeV 2 , is again quoted as 
an uncertainty. 

From ( |3.70| ) one can calculate the width for the decay D* — ► Dir. The prediction, 

T(D* + -> D\ + ) = I q I 3 = 32±5keV, (3.71) 

24:7im 2 Di , 
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Table 4: Theoretical estimates of the strong B*Bir and D*Dir couplings (from H2t\J). 



Reference g g B . B7l 9d*d-k T{D*+ -> D°tt+) (keV) 



29| 
29) Q 
l| a 
46l Q 



0.32 ± 0.02 29 ± 3 12.5 ± 1.0 32 ± 5 
28 ± 6 11 ± 2 
32 ± 6 

0.2 + 0.7 

0.39 ± 0.16 20 ± 4 9 ± 1 

0.21 ± 0.06 15 ± 4 7 ± 1 10 ± 3 

0.7 



0.75 1.0 
0.6 + 0.7 
0.4 + 0.7 
0.3 



64 



100 -i- 180 
61 78 



16.2 53.4 

19.5 ± 1.0 76 ± 7 

16.2 ± 0.3 53.3 ±2.0 

8.9 16 

8.2 13.8 

< 21 < 89 



a QCD sum rules in external axial field or soft pion limit. 

* including perturbative correction to the heavy meson decay constants. 

b Quark model + chiral HQET. 

c Chiral HQET with experimental constraints on D* decays. 

d Extended NJL model + chiral HQET . 

e Quark Model + scaling relation. 

' Relativistic quark model. 

9 Bag model. 

" SU(4) symmetry. 

1 Reggeon quark-gluon string model. 

k Experimental limits 
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lies well below the current experimental upper limit, 

T(D* + -> D\ + ) < 89 keV . 



(3.72) 



The latter is derived from the upper limit T tot (D* + ) < 131 keV and from the branching 
ratio BR(D* + — ► D°tt + ) = (68.3 ± 1.4) % |4"4]] . Predictions for other charge combinations 
are readily obtained from (|3.71|) and isospin relations analogous to ([3.671) . Accounting also 
for the differences in phase space, one expects 



r(£> H 



2.2T(D* + -> D + tt°) = 1.44r(D*°^D°7r 0> 



(3.73) 



In contrast to gD*Dn, the coupling constant Qb*b-k cannot be measured directly, since the 
corresponding decay B* — > Btt is kinematically forbidden. 
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Figure 11: The B — ► 7r /orm factor f + : direct sum rule prediction (solid), and single-pole 
approximation normalized by the sum rule estimate j or g^B-n (dashed). 

In Tab. 4, the numerical results discussed above are summarized and compared with 
other estimates. One observes significant differences. Some of the predictions are rather 
close to the experimental upper limit, some already violate the bound. It would be very 
interesting to have more precise data. 

As far as the asymptotic dependence of the VPn coupling on the heavy quark mass is 



concerned, the sum rule (|3.62|) suggests 



gvPn ~ m Q 



(3.74) 



in agreement with the expectation from HQET [^J, |48], [58| . The sum rules can also be used 
to investigate the 1/mg corrections. A simple quantitative estimate of the latter is obtained 
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p 2 [GeV z ] 



Figure 12: The D — > 7r /orm factor f + : direct sum rule prediction (solid), and single-pole 
approximation normalized by the sum rule estimate for go* Dn (dashed). 



by fitting the numerical predictions (|3.66| ) and ( [3.70D to the expression 



and the analogous one for go*D-K- The result is [p9| 



^ = 0.32 ±0.02, A = (0.7 ± 0.1) GeV 



(3.75) 



(3.76) 



In Tab. 4 , the above value of the reduced coupling constant g is compared to the estimates 
from other approaches. The 1/m.Q correction being about 15% for gs'B-K and 40% for go*D-K 
is non-negligible. Furthermore, in the heavy quark limit the ratio 

9B*B-KfB*\Jmr> , 

r = v (3.77) 

9d*d-kJd* a/ m B 

is expected to approach unity. Moreover, it has been shown []59| that r is subject to I/tuq 



corrections only in next-to-leading order. Indeed, the ratio r = 0.92 derived from the light- 
cone sum rules deviates surprisingly little from unity, in qualitative agreement with the 
HQET expectation. 

Finally, let us return to the original problem of calculating the heavy-to-light form factors 
at large momentum transfer. As pointed out in the beginning of this section, the VPir on- 
shell vertex is of great importance for the understanding of the form factor / + at large 
momentum transfer, since near the kinematic limit the V pole is expected to dominated the 
behaviour of / + . For B — > n, the single-pole approximation, 

fB*9B*Bir 



f + (p 2 



2m B *(l 



p 2 /m 2 



(3.78) 



B* j 
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is illustrated in Fig. 11 taking gn*B-K from (|3.66| ) and /b* from (|3.65|) . The extrapolation 
of the single-pole model to smaller p 2 matches quite well with the direct estimate from the 
light-cone sum rule (|3.41|) at intermediate momentum transfer p 2 = 15 to 20 GeV 2 . One 
thus has a consistent and complete theoretical description of / + . The extrapolation for the 
D — > 7r form factor using the analogous single-pole formula with go*Dn from (|3.70|) and fjj* 
from ( |3.69| ) is shown in Fig. 12. Also in this case we find the direct sum rule result and the 
pole model to match nicely at p 2 ~ 0.7 GeV 2 . However, here the matching point is a little 
bit too far away from the D* resonance for the extrapolation to be completely trustworthy. 



3.7 Decay spectra, integrated widths, and 

With the / + form factors for B and D mesons at hand, we are now in the position to 
discuss the differential distributions and integrated widths for the exclusive semileptonic 
decays B — > -nlv\ and D — > -nlv\ with I = e or /i. In these decays, the form factor f° being 
suppressed by mf plays a negligible role. 

For convenience, one may fit a parametrization of the form 



f + (P 2 ) = , ■_,/., 2 ,4/ 4 (3-79) 



/ + (0) 

1 — ap 2 /m 2 P + bp 4 /trip 

to the theoretical curves plotted in Fig. 11 and 12. For the B — > 7r form factor this fit yields 

a= 1.50, 6 = 0.52 . (3.80) 

Here, the known NLO corrections are included and, correspondingly, / + (0) has been fixed 
to the NLO value / + (0) = 0.27 given in ( |3.54| ). Fig. 13 shows the interpolation ( |3.79| ) in 
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Figure 14: Distribution of the momentum transfer squared in B — > -nlv\ (I = e,fi). The 
dashed curves indicate the theoretical uncertainty discussed in detail in the text. 



comparison with recent lattice results f\ The agreement is very encouraging. The analogous 
fit of (|3.79| ) for the D — > 7r form factor results in 

a = 1.16, b = 0.32 . (3.81) 

In this case, the NLO effects are not yet included whence the LO prediction / + (0) = 0.68 
given in ( |3.55| ) is used. 



The distribution in momentum transfer squared for B — > ixlvi is shown in Fig. 14. The 
band indicates the present theoretical uncertainty as discussed in detail in sect 3.5. To be 
precise, the uncertainties from the different sources are added in quadrature. Integrating 
this distribution one obtains the partial width 

F(B -> ixlvi) = (7.5 ± 2.5) \V ub \ 2 ps" 1 . (3.82) 

The corresponding distribution in the charged lepton energy is displayed in Fig. 15. 

Contrary to the semileptonic decays into e and /x, the decay B — ► irf v T is quite sensitive 
to the scalar form factor /°. As well-known [p9[| , the single-pole approximation used to 
extrapolate the sum rule result on / + to maximum p 2 cannot be applied to f°. One way 
to argue is that the scalar B ground state is about 500 MeV heavier than the pseudoscalar 
B and, therefore, lies too far above the kinematical endpoint of the B — > 7r transition in 
order to dominate the form factor. Nearby excited resonances and nonresonant states are 
expected to give comparable contributions. Thus, in order to demonstrate the sensitivity of 
the T-mode to f°, the light-cone prediction shown in Fig. 10 is linearly extrapolated from 



3 For a comprehensive review see J.M. Flynn and C.T. Sachraida, ref. 
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Figure 15: Distribution of the charged lepton energy in B — > nlui for I — e, \i (solid) and 
I = t (dashed). In the latter case the two curves correspond to the two extrapolations of f° 
shown in Fig. 16. 



the value at p 2 ~ 15 GeV 2 at which the sum rules (|3.41| ) and ( 3.42 ) still hold to the value 
at p 2 ~ m 2 B dictated by the soft pion limit |38|, |59fl : 

lim f{p 2 ) = If = 1.1 to 1.6 . (3.83) 

Here, we have used the conservative range of estimates fs = 150 to 210 MeV. This rough 
extrapolation of f° is illustrated in Fig. 16 pDf. Also shown are lattice data. They are 



systematically lower than the expectation. Note, however, that the sum rule estimate of 
f° still lacks NLO corrections. Certainly, within the theoretical uncertainties there is no 
contradiction. 

Fig. 15 shows the r-energy spectrum. As anticipated, a measurement of B — > 7ifu T has 
the potential to determine or at least constrain the scalar form factor. For the integrated 
width one predicts 

T(B nfu T ) = (6.1 ± 0.4)|Kfc| 2 PS" 1 • (3.84) 
Interesting is also the ratio of ( |3.84j ) and ( ^.82| ): 

T(B 7TfU T ) 



T(B — > 7rev e 



0.75 to 0.85 . (3.85) 



This expectation is independent of V u b, and less sensitive to uncertainties in the sum rule 
parameters than the widths themselves. The range quoted above corresponds to the two 
extreme extrapolations of f° considered. 
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Figure 16: TTie B ^ n form factor f°: direct sum rule estimate (solid) and linear extrapo- 
lations to the limit flff . 83j) (dashed). The lattice results are from fill 



Recently, the CLEO collaboration |2(J has reported the first observation of the semilep- 
tonic decays B — > 7r/z/ and -B — > p/z/ (Z = e,p). From the measured branching fraction, 
5^(5° -> 7T"Z+^) = (1.8 ± 0.4 ± 0.3 ± 0.2) • 10" 4 , and the world average of the B° lifetime 
[f44]|, r B o = 1.56 ± 0.06 ps, one derives 

T(B° -> tTZ+i/j) = (1.15 ± 0.35) • 10~ 4 ps" 1 , (3.86) 

where the errors have been added in quadrature. Comparison of (|3.86|) with Q3.82|) yields 

\V ub \ = 0.0039 ± 0.0006 ± 0.0006. (3.87) 

Here, the first (second) error corresponds to the current experimental (theoretical) uncer- 
tainty. 

A similar analysis has been performed for B — ► plv\ |^9[| . There are four independent 
B —>■ p form factors: 



(p(A)|m 7m (1 -7s)6| s > = ~*(™b + m p )A 1 (t)e(A) M + 



{e{\)p B )(pB +P P ) 



+ 



L4 3 (t) 



ms + m 



■{e{\)p B ){pB-P P \ 



(3.88) 



, 2Vft) £ ^ e(A) 



with t = (pb —Pp) 2 - Only V(i), ^i(t) and A 2 (i) contribute to the semileptonic decay into 
I = e or /i. The sum rule predictions on these form factors are plotted in Fig. 17, while Fig. 
18 shows the resulting distribution in momentum transfer squared. 
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Figure 17: JTie B —>■ p form factors: predictions from light-cone sum rules (solid) in com- 
parison to lattice results PdjJ . The dashed curves indicate the uncertainties in the sum rule 



results (from jjSj^J) 
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Figure 18: Distribution of the momentum transfer squared in B — > pev e . The dashed curves 
indicate the theoretical uncertainty (from WQD - 



From the integrated width [39 



T(B -> plvt) = (13.5 ± 4) \V ub \ 2 ps- 1 



and the CLEO result, 



T(B° -> /0"Z + Z4) = (1.60 ± 0.6) • 10~ 4 ps~ ] 



(3.89) 



(3.90) 



derived from the observed branching ratio BR(B° — ► p / + ^) = (2.5 ± 0.4lo;7 ± 0.5) ■ 10 4 
and the B° lifetime used in ( 3.86|) , one obtains 



ub 



0.0034 ±0.0006 ±0.0005. 



(3.91) 



Within errors, the values of \V u b\ extracted from the two exclusive measurements ( p. 86 ) and 
( |3.90|) are nicely consistent with each other, and also coincide with the inclusive determina- 
tion ofV ub |0. 



It is worth mentioning that the ratio 

r(£?o _ p-i+vj) 



1.8 ±0.7 



T(B° -> TT'l+Ul) 

following from ( |3.89| ) and ( p.82| ) is also consistent with the ratio 

BR(B° -> p -l + u e ) 



(3.92) 



1.4lg;4±0.3±0.4 



(3.93) 



observed by CLEO. However, the uncertainties on both sides are still too big to draw any 
firm conclusion. 
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The final application of the sum rule results on form factors which I want to consider 
here is the Cabibbo-suppressed decay D — > -neu e . The decay distribution ( |3.3| ) calculated 
from the parametrization (|3.79|) with (|3.81| ) for / + is plotted in Fig. 19. Since the electron 
mass can be neglected, the form factor / plays no role. The integrated width is predicted 
to be 

T(D -> ixlvi) = 0.16 \Ved\ 2 ps" 1 = 8.0 • 10~ 3 ps" 1 , (3.94) 

where |V^| = 0.224 ± 0.016 |44| has been used. This prediction should be compared with 
the experimental result, 



r(L> c 



(9-2i 2 2 ;l) • 10- 



ps 



(3.95) 



derived from the branching ratio BR(D° — ► n~e + u) = (3.8^;q) ■ 10~ 3 and the lifetime 
r D o = 0.415 ± 0.004 ps []44| . The theoretical uncertainty is estimated to be of the order 
of the experimental error. Despite of the preliminary agreement, one cannot be satisfied 
with the present status. In order to really test the sum rule approach both more accurate 
calculations and more precise data are needed. 
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Figure 19: Distributions of the momentum transfer squared in D — » 7reu e . 



4 Nonleptonic two-body decays 
4.1 Effective hamiltonian 

As illustrated in Tab. 2, exclusive nonleptonic decays are the theoretically most complicated 
processes discussed in these lectures. This class of decays is strongly influenced by hard and 
soft interactions of quarks and gluons, by initial bound state and nonspectator effects, as well 
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as by hadronization and final state interactions. Although over the years one has developed 
a qualitative understanding and even an amazingly consistent phenomenological description 
of two-body decays [], some features are still lacking a quantitative theoretical explanation. 

Up to now, only the hard-gluon effects can be systematically taken into account in 
the framework of renormalization group improved QCD perturbation theory. As already 
pointed out, the essential tool is OPE which allows to separate the dynamics at short and 
long distances. The result is an effective weak hamiltonian given by a sum of renormalized 
local operators with scale-dependent Wilson coefficients as sketched in ( p..2|) . Specifically 
for nonleptonic B decays, one has 

Q 

H NL =-^ V qb [ Cl (fi)Oi q (n) + c 2 (n)0 2q (ii)] + h.c. , (4.1) 

V " q=c,u 

where 

lq = (d'T^u + rs'^c)(qT p b), 2q = {qT»u){d'T p b) + (qVc) (s'T p b) (4.2) 

with r p = 7 p (l — 75 ). Similarly, the effective hamiltonian relevant for nonleptonic D decays 
is given by 

G 

H NL = ^=Y. V «} [cMOiq(v) + c 2 (n)0 2q (fi)} + h.c. , (4.3) 

V 2 q=s4 

with 

lq = {d'TPu){cT p q), 2q = (cn)(d'T p q) . (4.4) 

In the above, d! = V u dd + V us s and s' = V c dd + V cs s are weak eigenstates. In addition to the 
four-quark operators 0\ and 2 , the operator-product expansion includes so-called penguin 
operators H63| . These are mainly relevant for rare nonleptonic decays such as B — > Kit or 
radiative decays such as B —>■ K*j, but they play no role in the following discussion. 

The Wilson coefficients Cj(/i) contain the effects from hard gluon and quark exchanges 
with virtualities larger than \i. In perturbation theory, these interactions generate large 
logarithmic terms, a™Zn n! ^, which must be summed up to all orders. This is achieved by 
solving the renormalization group equation 

(4.5) 
where 

a s (i) / a s \ 2 (2) 
7± = 47 7± + UJ 7± + (4 ' 6) 

are the anomalous dimensions of the operators 0± = |(Oi±0 2 ), and ci j2 = |(c + ±c_). With 
the boundary conditions c±(m\y) = 1 following from the unrenormalized weak hamiltonian 
indicated in ( |1.1| ) the solution is given by 

c ± (/,) = ( ^L) ^ (l + ^ "-M - °>(™*) ) . (4 .7) 
\a s (m w ) I \ 4vr I 




3 For a recent comprehensive review see M. Neubert and B. Stech, ref. [B2| 
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The first factor is the LO result |M , while the second bracket incorporates the NLO cor- 
rections |65|. The LO coefficients are determined by the one-loop anomalous dimensions, 



7? 



+8, 7 



and by the one-loop coefficient of the QCD /3-function, 



6=11- -n f 



(4.8) 



(4.9) 



The NLO term R± is given in |)5] 



The relevant physical scale is expected to be of the order of the heavy quark mass. For 
ji = m b = 4.8 GeV one obtains p2|, numerically, 



while for /i 



Ci(m 6 ) = 1.108, c 2 (m b ) 

ci{m b ) = 1.132, c 2 {m b ) = 
1.4 GeV one gets 

Ci(m c ) = 1.263, c 2 (m c ) 

Ci{m c ) = 1.351, c 2 (m b ) = 

Here, the running coupling is normalized to a s (mz) = 0.118. As anticipated, the effects are 
smaller in b decays than in c decays. Furthermore, the NLO effects re-enforce the LO trend. 



-0.249 (LO) , 
-0.286 (NLO) . 

-0.513 (LO), 
-0.631 (NLO) 



(4.10) 
(4.11) 

(4.12) 
(4.13) 



4.2 Decay amplitude and factorization 

Given the effective hamiltonian, the amplitude for a two-body decay is calculated from the 
corresponding hadronic matrix elements of the local four-quark operators: 

A{P^h 1 h 2 ) = {h 1 h 2 \H NL \P) 

= %E^E^(/i)(M 2 |0 4 »|P). (4.14) 

V2 q i 

Since these matrix elements contain soft QCD interactions characterized by virtualities 
smaller than /i, the problem of calculating them is extremely difficult and still far from a 
satisfactory solution. It may not be completely needless to stress that the scale dependence 
of the matrix elements must cancel in ( f4 . 1 4| ) against the scale dependence of the Wilson 
coefficients since the decay amplitude is an observable. In practice, there is always some 
residual scale dependence because of the truncation of the perturbative series. In terms of 
the weak amplitude ( [4.1 4|) the partial width for two-body decays is given by 



T(P -> hM = t\ (ml + m? - m 2 2 ) 2 - Amlml \A\ 2 . (4.15) 

167rm; ' 



up 



The present knowledge and understanding of nonleptonic B decays has recently been 
reviewed in a very comprehensive article by M. Neubert and B. Stech ||62|| . For D decays 



one may consult the classic papers by M. Bauer, B. Stech and M. Wirbel ]4T|. Here, I will 
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restrict the discussion to a particular example which brings the main theoretical difficulties 
to light, namely B — > J/ipK. This decay mode also plays an outstanding role in I?-physics 
at hadron colliders and future 5-factories. 

After Fierz-transformation of the operator 0\ c in (43) the piece of the effective hamil- 
tonian relevant for B — > J/ipK can be written in the form 

H NL = -^=V cb V; s {(c 2 + |)0 2 + 2 Cl 2 } (4.16) 

with 

2 = (cT p c)(sT p b), 6 2 = (cT^c)(sT p ^-b) , (4.17) 

A a being the usual SU(3)-co\ovlt matrices. The colour indices and the /i-dependence are 
suppressed for simplicity. 

In a radical first approximation, one may factorize the matrix elements of the operators 
2 and 2 into products of matrix elements of the currents that compose these operators. 
Long distance strong interactions between quarks belonging to different currents are thereby 
completely neglected. Furthermore, the matrix element of the octet operator 2 vanishes 
because of colour conservation, while the factorized matrix element of the singlet operator 
2 is approximated by 

(J/iPK | 2 (/i) | B) = (J/tfj | cT p c \0)(K\ sT p b | B) = 2f i) f+_ K m^e >l) ■ q) . (4.18) 

Obviously, is the J/ip decay constant, Jb^k * s ^ ne B ^ K form factor at the momentum 
transfer p 2 = m^, denotes the J/if> polarization vector, and q the K four-momentum. 
The leptonic width T(J/ip^ = 5.26 ± 0.37 keV implies 

= 405 MeV , (4.19) 

while a sum rule estimate |28j similar to the one for the B — > ir form factor in sect. 3.5 
yields 

f+^ K = 0.55 ± 0.05 . (4.20) 

At this point one encounters a first principal problem: since the matrix elements of 
the quark currents in (|4.18|) are scale- independent, the /i-dependence of the short-distance 
coefficients cannot be compensated. Hence, the approximation 

(J/xPK | H NL | B) = V2G F V cb V c * s (cM + ^j^j Uf^K™^ ■ o) (4-21) 

can at best be valid at a particular scale fi which could be called the factorization scale. 
The conventional assumption is fi = 0(m&). 

Using the next-to-leading order coefficients Ci j2 (/i) in the HV scheme with = 225 
MeV from |)6| and taking /x = m h ~ 5 GeV, one has 

c 2 (/i) + = 0.156 . (4.22) 

Together with (|4.19 ) and ( 4.20|) , this yields the branching ratio 

BR(B -> J/ipK) = 0.025% , (4.23) 
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a value which is considerably smaller than the measured branching ratios f44] , p7| 

BR(B- -> J/ipK') = (0.101 ± 0.014)% , (4.24) 

BR(B° -> J/^K°) = (0.075 ± 0.021)% . (4.25) 

The discrepancy between experiment and expectation is worsened by a factor 3 if the coef- 
ficients ( [4.111) are used yielding c 2 + ci/3 = 0.091 instead of Q4.22Q . This provides a good 
example for the scale uncertainty due to factorization. 

4.3 Nonfactorizable contributions and effective coefficients 

The quantitative failure and the scale problem pointed out above imply that naive factor- 
ization of matrix elements does not work. Clearly, factorization has to be accompanied by 
a reinterpretation of the Wilson coefficients. For B — > J/ipK, the short-distance coefficient 
c 2 (/x) + Ci(/i)/3 is substituted by an effective scale-independent coefficient a 2 . Phenomeno- 
logically f4"I|| , a 2 is treated as a free parameter to be determined from experiment. From 
flPID and ( pL2p one finds 



jaf^l = 0.31 ±0.02 , (4.26) 

where the quoted error is purely experimental. The sign of remains undetermined. It 
is an interesting observation that the above value is similar to the coefficient a 2 found from 
so-called class III decays p2fl . This suggests, but does not prove a universal nature of a 2 . 



The drastic difference between the short-distance coefficient (|4.22|) and the effective 



coefficient ( |4.26| ) points at the existence of sizeable nonfactorizable contributions. The 



latter are also needed to cancel or at least soften the strong //-dependence of the factorized 
amplitude ( |4.21|) . Further investigation shows that the dominant nonfactorizable effects 



should arise from the matrix element of the operator 2 . Writing the latter in the convenient 
parametrization 

(J/i>K | 6 2 (ji) | B) = 2Uf B ^ K {p)m^ ■ q) , (4.27) 
and adding it to ( [4.18| ), one formally reproduces ( [4.2 1| ), 

(J/i/>K | H NL | B) = V2G V cb V; s a^ K Uf^ K m4e^ ■ q) , (4.28) 

with the effective coefficient |68|] 

ag+ K = cM + °-^ + 2 Cl (/,)%^ . (4.29) 

^ JB^K 

In Fig. 20, the partial width for B — > J/ipK is shown as a function of the parameter fstpK 
associated with (0 2 ). Note that /b^k = corresponds to naive factorization, while the 
fitted value of a 2 ^ K given in ( [4.26| ) implies 

fwiciti = m b ) = +0.04 or - 0.12 . (4.30) 

For a lower scale the value of Jb^k is shifted to slightly more positive values, e.g., 

hMv = \ m b) = +0-06 or - 0.09 . (4.31) 

We see that a nonfactorizable amplitude of 10 to 20 % of the size of the factorizable one is 
sufficient to reconcile expectation with experiment. 
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T(B -> J/ipK) [10-V 
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Figure 20: The partial width for B — > J/ipK as a function of fsipK parameterizing the 
nonfactorizable contribution to the decay amplitude. The horizontal dotted lines represent 
the average of the experimental widths given in the text. The solid (dashed) curves show 
the theoretical expectation for // = rrib/2 (fi = nib). The arrows indicate the QCD sum rule 
estimate. 



4.4 Sum rule estimate of the nonfactorizable amplitude 

The calculation of a^ K and of the analogous coefficients for other two-body decays |4l| 



from first principles is one of the outstanding problems in weak decays. As a first step 
in this direction, we have undertaken a rough estimate of fs^K using again QCD sum 
rule methods Q Following the idea worked out for charm decays in [BS], we consider the 
four-point correlation function 

n^(p, q) = J d'x d 4 y d'z e ^ + ^(0 | T{j«(x)ji(y)d 2 (z)ji(0)} | 0) , (4.32) 

where jjf = bij^u, jf = c^ u c, and = u^/y§s are the generating currents of the mesons 
participating in the decay B — > J / ipK, and p + q, p and q are the respective four- momenta. 
O2 is the octet operator given in (|4.17|) which is suspected to give rise to the dominant 
effect. 

Similarly as in the case of the two-point correlation function (|2.4|) studied in sect. 2.2, 
one writes a dispersion relation for (4.32) in terms of intermediate hadronic states in the B, 



7 This work was done in collaboration with B. Lampe and A. Khodjamirian. Results have been reported, 
e.g., in f38|. 
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J/i/j, and K channel. The ground state contribution to (|4.32 



n ( v _ . (0 | 3% | K)(0 | jf | J/jj)(J/jjK | Q 2 | B)(B \ jf | 0) 
VlP ' gj "' «-g 2 )K-p 2 )(m|-(p + g) 2 ) + ' ^ 66) 

contains the desired matrix element ( |4.27|) . In addition, one has contributions from excited 
resonances and continuum states denoted above by ellipses which lead to a very compli- 
cated singularity structure. Furthermore, in the euclidean momentum region Q 2 = —q 2 > 
0(1 GeV 2 ), p 2 < 0, (p + q) 2 < one can expand the correlation function Q4.32|) in terms of 
local operators: 

n^(p,g) = Y,C?ip,q,tiW)) ■ (4-34) 



For the operators with dimension d < 6 given in ( |2.13| ), one can calculate the necessary 



coefficients C^ u (p,q, /a) from the diagrams shown in Fig. 21. From (|4.33| ) and (|4.34j ) it is 
possible, at least in principle, to derive a sum rule for the matrix element (J/ipK | 2 | B). 

However, there are two complications that are not present in the two point sum rules 
discussed in sect. 2.2. One problem is the presence of a light continuum in the B channel 
below the pole of the ground state 5-meson. This contribution to ( }4.33|) can be associated 



with processes of the type B — > "D*D S " — > J/ipK, where an intermediate state carrying 
D*D S quantum numbers rescatters into the final J/ipK state. Formally, in ( 4.32| ) this 



intermediate state is created from the vacuum by the combined action of the operator 
product 2 j^. As a reasonable solution we suggest to cancel this unwanted contribution 
against the terms in ( f4.34|) which have the quark content ccsq and which develop a nonzero 
imaginary part at (p + q) 2 > 4m 2 . In the approximation considered this is the case for the 
four-quark condensate contribution represented by Fig. 21c. 

The second problem concerns the continuum subtraction in the remainder. A detailed 
examination of (|4.32|) indicates that the higher charmonium resonances contribute with 
alternating signs. Therefore, the usual subtraction procedure in which the dispersion integral 
over the excited and continuum states is approximated by the perturbative counterpart is not 
reliable here. In order to proceed we employ explicit, although rough models for the hadronic 
spectral functions. Since the number of additional parameters has to be manageable, only 
the first excited resonances are included in each channel besides the B, J/ip and K ground 
states. In total, in this approximation the correlator (f4.33|) contains three free parameters 



in addition to fsipK- Finally, one takes the Borel transform of ( }4.33|) and Q4.34 ) in the 
5-meson channel and moments in the charmonium channel. In the i^-meson channel, q 2 is 
kept spacelike. 

A fit of the hadronic representation ( 4.33 ) to the OPE result ( [4.34 ) for various values of 
the Borel mass M and q 2 , and for different moments yields 

f Bi , K = -(0.045 to 0.075). (4.35) 



The implicit scale of this estimate is given by M ~ ymjs ~ — \rrib — 2.4 GeV. Substi- 
tuting (|4.35|) in (|4.29|) , and evaluating the short-distance coefficients Ci^fj) also at // = M 



one gets 

a B^K = _ Q 2g + 3g _ ( ig fo 31 ^ = _( .io to 0.22) . (4.36) 

Here, the three terms in the first relation refer to the three terms in (|4.29| ) in the same 
order. Interestingly, the sum rule approach seems to favour the negative solution for fsipK- 
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Although in comparison with Q4.31 ) the above estimate falls somewhat short, it narrows the 
gap between theory and experiment considerably as can be seen from Fig. 20. 





Figure 21: Diagrams determining the Wilson coefficients of the OPE of the correlation 
function 11^. The symbols are as in Fig. 1. 

Several comments are in order. Firstly, the nonfactorizable matrix element ( f4.27p is 
rather small as compared to the factorizable one given in ( |4.18| ), numerically, 



\Ib^k/ f, 



0.1 . 



(4.37) 



Nevertheless, the impact on af ^ K is strong because of the large coefficient as compared to 
the coefficient of the factorized amplitude, 



|2ci/(c 2 + ci/3)| ~20 to 30 



(4.38) 



Secondly, the factorizable amplitude proportional to Ci/3 and the nonfactorizable one pro- 
portional to fs^K are opposite in sign and hence tend to cancel. In fact, if \ fs^K\ is taken 
at the upper end of the estimated range, the cancellation is almost complete resulting in a 
considerable enhancement of the branching ratio ( ^.23| ) . Note that both terms are nonlead- 
ing in 1/N C . This is exactly the scenario anticipated by the 1/N C - rule for D decays |70|] 
which is generally consistent with observation |4j] and has found theoretical support by the 
sum rule analysis of [B9|. In |7T| , 72 1, a similar trend was claimed for B — > Dir. Thirdly, 
our estimate yields a negative overall sign for a^ K in contradiction to a description of 



45 



the data on two-body B decays based on factorized decay amplitudes with two universal 
coefficients a\ and a 2 |62]| . It should be stressed, however, that in this fit the positive sign 
of a 2 actually results from the channels B~ — > D°tt~, D°p~ , D*°tt~ , and D*°p~ , and is then 
assigned also to the J/ipK channel. This assignment may not be correct. Certainly, the 
sum rule approach described above provides no justification for such an assumption. On the 
contrary, diagrams of the kind shown in Fig. 21 suggest some channel-dependence of the 
nonfactorizable matrix elements. For class II processes the channel-dependence is enhanced 
by the large coefficient 2ci in a 2 (see ( |4.29| )). The opposite is the case for class I processes 
involving 

a 1 =c 1 (/i) + ^ + 2c 2 -^ : . (4.39) 

Here, the nonfactorizable contributions are damped by the small coefficient c 2 , while the 
factorized term has a large coefficient. Therefore, a\ is indeed expected to be universal to 
a good approximation, whereas a 2 should exhibit some channel-dependence, in particular 
when comparing decays with very different final states such Dtt and J/ipK. 

As a last remark, from the sum rule point of view one does not expect a simple relation 
between B and D decays. At least, the OPE of the four-point correlation functions (such as 
fl4.32| )) is significantly different, e.g., for the channels B — > Dti, B — > J/ipK, and D — > Ktt 
as can be inferred from Fig. 21. Hence, the arguments presented in |)2| in favour of a 
change of sign in ct 2 when going from D to B may oversimplify the dynamics of nonleptonic 
decays. 



5 Conclusion 



These lectures have been devoted to the theory of exclusive decays of B and D mesons. I 
have discussed examples of leptonic, semileptonic and nonleptonic decays, focussing on the 
problem of calculating hadronic matrix elements of weak currents and four-quark operators. 
The aim has been to demonstrate the usefulness of QCD sum rule techniques. I have 
described the derivation of sum rules for decay constants, form factors, and two-body decay 
amplitudes. In addition, I have outlined a sum rule for couplings between heavy and light 
mesons. Numerical predictions have been presented for the decay constants fo, /d 3 , Ib, and 
f Bs , for the heavy-to-light formfactors f + (p 2 ), f~(p 2 ), and f°(p 2 ), for the strong couplings 
9b*b-k and gD*Dn, and for the nonfactorizable matrix element {J/ipK | 2 (/x) | B). 

Using the sum rule results I have presented predictions on decay distributions and in- 
tegrated widths for B — > ttIvi, B — > plui, and D — > ixlvi. Comparison with the CLEO 
measurements |20| of exclusive semileptonic B decays yields values for V u b in good agree- 
ment with each other and with the determination from inclusive data Agreement 



between expectation and measurement is also found for the Cabibbo-suppressed semilep- 
tonic D° decay |44j]. Furthermore, a sum rule estimate is given for D* — > Dir. The present 
experimental upper limit is still about three times larger than the expected width. Fi- 
nally, including the sum rule estimate of the nonfactorizable contribution to the amplitude 
for B — > J/ipK a theoretical estimate of the effective coefficient a 2 is obtained which is 
rather close to the value extracted from the experimental branching ratio pi |. However, 
the sign of a 2 predicted by the sum rule analysis is opposite to the sign determined from 
data, if channel-independent, universal coefficients ai j2 are assumed fl62|. I have explained 
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why in the sum rule approach universality can be expected for ai, but not for 0,2- From the 
sum rule point of view the apparent universality of a\ and a<± constitutes a major puzzle 
which needs further clarification. 

Finally, I have addressed the present theoretical uncertainties in the sum rule calcula- 
tions, and the prospects for reducing the uncertainties. On the theoretical side, one has 
to improve the light-cone wave functions of light mesons, and calculate the higher-order 
perturbative corrections to the sum rules. On the experimental side, new and more precise 
measurements at hadron facilities and at future B and tau-charm factories should allow to 
tightly constrain the input parameters and to test the reliability of the sum rule approach. 
It appears conceivable to ultimately decrease the uncertainties from presently 20 to 30 % 
by a factor of 2. Encouragement comes from the preliminary agreement of lattice and sum 
rule calculations on decay constants and form factors. It would certainly be very fruitful 
to combine the flexibility of the sum rule method with the rigorous nature of the lattice 
approach. 

Acknowledgements 

I want to thank Ikarus Bigi and Luigi Moroni for having invitated me to the School 
"Enrico Fermi", and the staff of the school for having made my stay so enjoyable. I am 
particularly grateful to Alexander Khodjamirian for his continuous collaboration and for 
his great help with the present manuscript. Furthermore, I thank V.M. Belyaev, V.M. 
Braun, B. Lampe, Ch. Winhart, S. Weinzierl, and O. Yakovlev for collaboration on various 
subjects covered in these lectures and for useful discussions. This work was supported by the 
Bundesministerium fiir Bildung, Wissenschaft, Forschung und Technologie, Bonn, Germany, 
Contract 05 7WZ91P (0). 

References 

URALTSEV N.G., this volume. 



[2 

[3 

[4 

[5 
[6 

[7; 



SHIFMAN M.A., VAINSTHTEIN A.I. and ZAKHAROV V.I., Nucl. Phys. B, 147 
(1979) 385; 448. 

REINDERS L.J., RUBINSTEIN H.L. and YASAKI S., Phys. Rep., 127 (1985) I. 

Vacuum Structure and QCD Sum Rules, edited by M.A. SHIFMAN (North-Holland, 
Amsterdam) 1992. 

KHODJAMIRIAN A. and RUCKL R., to appear in Heavy Flavours, 2nd edition, edited 
by A.J. BURAS and M. LINDNER (World Scientific, Singapore) 1998, [hep-ph/9801443l 



RICHMAN J.D., in Proc. of 28th ICHEP (Warsaw), edited by Z. AJDUK and A.K. 
WROBLEWSKI (World Scientific, Singapore) 1996, 143. 

ADLER J. et al. (MARK III Collab.), Phys. Rev. Lett, 60 (1988) 1375. 

BABAR Technical Design Report, BOUTIGNY D. et al. (BABAR Coll.), SLAC-R-95- 
457. 



47 



[9] NOVIKOV V.A., OKUN L.B., SHIFMAN M.A., VAINSHTEIN A.I., VOLOSHIN M.B. 
and ZAKHAROV V.I., Phys. Rev. Lett, 38 (1977) 626. 

[10] BROADHURST D.J., Phys. Lett. B, 101 (1981) 423; 

BROADHURST D.J. and GENERALIS S.C., preprint OUT-4102-8/R , 1982 (unpub- 
lished) . 

[11] ALIEV T.M. and ELETSKY V.L., Sov. J. Nucl. Phys., 38 (1983) 936. 

[12] VOLOSHIN M.B., Int. J. Mod. Phys. A, 10 (1995) 2865; 
JAMIN M. and PICH A., Nucl. Phys. B, 507 (1997) 334; 
KUHN J.H., PENIN A. A. and PIVOVAROV A. A., |hep-ph/9801356 . 

[13] GASSER J. and LEUTWYLER H., Phys. Rep., 87 (1982) 77. 

[14] BELYAEV V.M. and IOFFE B.L., Sov. Phys. JETP, 56 (1982) 493. 

[15] DOMINGUEZ C.A., Talk at 3rd Workshop on the Tan-Charm Factory (Marbella, 
Spam), 1993, [hep-ph/ 9309260) . 

[16] FLYNN J.M., in Proc. of 28th ICHEP (Warsaw), edited by Z. AJDUK and A.K. 
WROBLEWSKI (World Scientific, Singapore) 1996, p. 335, |hep-lat/96110T6 . 



[17] ALLTON C.R., CONTI L., CRISAFULLI M., GIUSTI L., MARTINELLI G. and RA- 
PUANO F., Phys. Lett. B, 405 (1997) 133. 

[18] WITTIG H., Int. J. Mod. Phys. A, 12 (1997) 4477. 

[19] DOMINGUEZ C.A. and PAVER N., Phys. Lett. B, 197 (1987) 423; 
NARISON S., Phys. Lett. B, 198 (1987) 104; 
REINDERS L.J., Phys. Rev. D, 38 (1988) 423. 

[20] ALEXANDER J.P et al. (CLEO Collab.), Phys. Rev. Lett, 77 (1996) 5000. 

[21] AITALA E.M. et al. (Fermilab E791 Collab.), Phys. Rev. Lett, 80 (1998) 1393; 

FRABETTI PL. et al. (Fermilab E687 Collab.), Phys. Lett. B, 364 (1995) 127; 382 
(1996) 312; 391 (1997) 235; 

BARTELT J. et al. (CLEO Collab.), Phys. Lett. B, 405 (1997) 373; 
BEAN A. et al. (CLEO Collab.), Phys. Lett. B, 317 (1993) 647. 

[22] CHERNYAK V.L. and ZHITNITSKY A.R., JETP Lett, 25 (1977) 510; Yad. Fiz., 31 
(1980) 1053; 

EFREMOV A.V. and RADYUSHKIN A.V., Phys. Lett. B, 94 (1980) 245; Teor. Mat. 
Fiz., 42 (1980) 147; 

LEPAGE G.P and BRODSKY S.J., Phys. Lett. B, 87 (1979) 359; Phys. Rev. D, 22 
(1980) 2157. 

[23] BRODSKY S.J. and LEPAGE G.P., in Perturbative Quantum Chromodynamics, edited 
by A.H. MUELLER (World Scientific, Singapore), 1989, 93. 

[24] CHERNYAK V.L. and ZHITNITSKY A.R., Phys. Rep., 112 (1984) 173. 



48 



[25] BALITSKY I.L, BRAUN V.M. and KOLESNICHENKO A.V., Sov. J. Nucl. Phys., 44 
(1986) 1028; Nucl. Phys. B, 312 (1989) 509. 

[26] BRAUN V.M. and FILYANOV I.B., Z. Phys. C, 44 (1989) 157. 

[27] CHERNYAK V.L. and ZHITNITSKY I.R., Nucl. Phys. B, 345 (1990) 137. 

[28] BELYAEV V.M., KHODJAMIRIAN A. and RUCKL R., Z. Phys. C, 60 (1993) 349. 

[29] BELYAEV V.M., BRAUN V.M., KHODJAMIRIAN A. and RUCKL R., Phys. Rev. 
D, 51 (1995) 6177. 

[30] KHODJAMIRIAN A., RUCKL R. and WINHART C, [hep-ph/ 9802412 . 

[31] BALITSKY I.L and BRAUN V.M., Nucl. Phys. B, 311 (1988) 541. 

[32] GORSKY A.S., Sov. J. Nucl. Phys., 41 (1985) 1008; ibid., 45 (1987) 512; ibid., 50 
(1989) 498 . 

[33] BRAUN V.M. and FILYANOV LB., Z. Phys. C, 48 (1990) 239. 

[34] KHODJAMIRIAN A., RUCKL R., WEINZIERL S. and YAKOVLEV O., Phys. Lett. 
B, 410 (1997) 275. 

[35] BAGAN E., BALL P. and BRAUN V.M., Phys. Lett. B, 417 (1998) 154. 

[36] DITTES F.M. and RADYUSHKIN A.V., Phys. Lett. B, 134 (1984) 359; 
SARMADI M.H., Phys. Lett. B, 143 (1984) 471; 

MIKHAILOV S.V. and RADYUSHKIN A.V., Nucl. Phys. B, 254 (1985) 89. 

[37] ISGUR N. and WISE M.B., Phys. Rev. D, 41 (1990) 151; ibid., 42 (1990) 2388. 

[38] VOLOSHIN M.B., Sov. J. Nucl. Phys., 50 (1989) 105. 

[39] BALL P. and BRAUN V.M., Phys. Rev. D, 55 (1997) 5561. 

[40] BALL P., BRAUN V.M. and DOSCH H.G., Phys. Lett. B, 273 (1991) 316. 

[41] WIRBEL M., STECH B. and BAUER M., Z. Phys. C, 29 (1985) 637; 
BAUER M., STECH B. and WIRBEL M., Z. Phys. C, 34 (1987) 103. 

[42] COLANGELO P. and SANTORELLI P., Phys. Lett. B, 327 (1994) 123. 

[43] ALI A., BRAUN V.M. and SIMMA H., Z. Phys. C, 63 (1994) 437. 

[44] Particle Data Group, Phys. Rev. D, 54 (1996) 1. 

[45] OVCHINNIKOV A.A., Sov. J. Nucl. Phys., 50 (1989) 519. 

[46] GROZIN A.G. and YAKOVLEV O.L, preprint BUDKERINP-94-3, 
|hep-ph/940 1267| (unpublished). 



49 



[47] COLANGELO P., NARDULLI G., DEANDREA A., DI BARTOLOMEO N., GATTO 
R. and FERUGLIO F., Phys. Lett. B, 339 (1994) 151. 

NUSSINOV S. and WETZEL W., Phys. Rev. D, 36 (1987) 130. 

YAN T.-M., CHENG H.-Y., CHEUNG C.-Y., LIN G.-L., LIN Y.C. and YU H.-L., 
Phys. Rev. D, 46 (1992) 1148, Erratum, ibid., 55 (1997) 5851. 

CHO P. and GEORGI H., Phys. Lett. B, 296 (1992) 408. 

AMUNDSON J.F., BOYD C.G., JENKINS E., LUKE M., MANOHAR A.V., ROSNER 
J.L., SAVAGE M.J. and WISE M.B., Phys. Lett. B, 296 (1992) 415. 

BARDEEN W.A. and HILL C.T., Phys. Rev. D, 49 (1994) 409. 

EICHTEN E., GOTTFRIED K., KINOSHITA T., LANE K.D. and YAN T.M., Phys. 
Rev. D, 21 (1980) 203. 

O'DONNELL P.J. and XU Q.P, Phys. Lett. B, 336 (1994) 113. 
MILLER G.A. and SINGER P., Phys. Rev. D, 37 (1988) 2564. 
THEWS R.L. and KAMAL A.N., Phys. Rev. D, 32 (1985) 810. 
KAIDALOV A.B. and NOGTEVA A.V., Sov. J. Nucl. Phys., 47 (1988) 321. 
WISE M.B., Phys. Rev. D, 45 (1992) R2188. 

BURDMAN C, LIGETI Z., NEUBERT M. and NIR Y., Phys. Rev. D, 49 (1994) 2331. 
FLYNN J.M. and SACHRAIDA C.T., |hep-lat/9710057| . 

Bid I., SHIFMAN M. and URALTSEV N., Ann. Rev. Nucl. Part. Sci., 47 (1997) 591; 
NEUBERT M., |hep-ph/9702375 . 



NEUBERT M. and STECH B., frep-ph/9705252 



SHIFMAN M.A., VAINSHTEIN A.I. and ZAKHAROV V.I., Nucl. Phys. B, 120 (1977) 
316. 

ALTARELLI G. and MAIANI L., Phys. Lett. B, 52 (1974) 351; 
GAILLARD M.K. and LEE B.W., Phys. Rev. Lett, 33 (1974) 108. 

ALTARELLI G., CURCI C, MARTINELLI G. and PETRARCA S., Phys. Lett. B, 99 

(1981) 141; Nucl. Phys. B, 187 (1981) 461; 

BURAS A.J. and WEISZ PH., Nucl. Phys. B, 333 (1990) 66. 

BURAS A.J., Nucl. Phys. B, 434 (1995) 606. 

ALAM M.S. et al. (CLEO Collab.), Phys. Rev. D, 50 (1994) 43. 

KHODJAMIRIAN A. and RUCKL R., Nucl. Instrum. Methods A, 368 (1995) 28. 



50 



[69] BLOK B.Yu. and SHIFMAN M.A., Sov. J. Nucl. Phys., 45 (1987) 135, 301, 522. 

[70] BURAS A.J., GERARD J.-M. and RUCKL R., Nucl. Phys. B, 268 (1986) 16. 

[71] BLOK B. and SHIFMAN M., Nucl. Phys. B, 389 (1993) 534. 

[72] HALPERIN L, Phys. Lett. B, 349 (1995) 548. 



51 



